Warping effects in the band and angular-momentum structures of the topological insulator Bi2Te3 
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We performed angle-resolved photoemission (ARPES) experiments on Bi2Te3 with circularly polarized light. 
ARPES data show very strong circular dichroism, indicating the existence of orbital angular momentum (0AM). 
Moreover, the alignment of 0AM is found to have a strong binding energy dependence. Such energy dependence 
comes from a relatively strong band warping effect in Bi2Te3 compared to Bi2Se3. 0AM close to the Dirac 
point has an ideal chiral structure (sin 0) without an out-of-plane component. The warping effect comes in as the 
binding energy decreases, and circular dichroism along a constant energy contour can no longer be explained 
by a simple sin 6 function but requires a sin W term. When the warping effect becomes even stronger near the 
Fermi energy, circular dichroism gains an additional sin 60 term. Such behavior is found to be compatible with 
the theoretically predicted 0AM structure. 
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Topological insulators (TIs) are a new class of insulators. 
They are distinguished from ordinary insulators by the Z2 
topological number. The most important aspect of such a 
topological state is the existence of topological metallic (TM) 
states at the boundary, that is, on the surface. These metallic 
states have Dirac dispersions and are protected against time- 
reversal symmetry-conserving perturbations. In addition, the 
states on the opposite sides of the Dirac cone are time-reversal 
pairs, resulting in suppression of backscattering. The time- 
reversal pairs are normally realized in the form of opposite 
electron spin states due to the Rashba effect from spin-orbit 
coupling, forming chiral spin states along a constant energy 
contour. 

In an ideal case, the Dirac dispersion of TM states 
is linear and has circular constant energy contours in the 
momentum space. ^ In such a case, backscattering of the 
surface electrons by nonmagnetic impurities should be com- 
pletely suppressed.^'^ However, angle-resolved photoemission 
(ARPES) data from Bi2Se3 and Bi2Te3 show hexagonal 
warping of the surface state Fermi surfaces. ^^'^^ Subsequent 
theoretical studies reveal that a strong warping effect can 
modify the spin chirality and brings in an out-of-plane spin 
polarization component. ^^'^^ Considering the fact that the 
spin chirality on the Fermi surface can affect the transport 
properties, it is desirable to have systematic experimental 
investigations of the warping effect on the spin states. 
Unfortunately, detection of the spin state by spin-resolved 
photoemission is not an easy task due to its low efficiency. 

Recent ARPES data that utilized circularly polarized light 
(CPL) may provide a way to circumvent the problem. The 
data from Bi2Se3 show a large circular dichroism (CD), 
which stems from the existence of orbital angular momentum 
(0AM). It was found that the QAM is also locked to the 
electron momentum and its direction is antiparallel to the 
spin angular momentum (SAM) direction. Therefore QAM 
also forms a chiral structure and is an essential ingredient 
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previously missing in the description of the surface Dirac 
states. Such observation provides an opportunity to find the 
spin direction (by measuring the QAM direction) without spin 
detectors, which have very low efficiency. One may thus be 
able to investigate the warping effect on the spin structure in 
Bi2Te3 by using dichroic ARPES. 

We report on our attempt to investigate the warping effect 
on the spin structures of the surface Dirac states of Bi2Te3. Our 
main goal was to obtain the systematic change in the 0AM 
(hence the SAM, too) structure in the momentum space at 
a constant binding energy as a function of the energy from 
the Dirac point. To achieve the goal, we performed ARPES 
experiments with circularly polarized light and obtained 
dichroism data. We observe the 0AM direction from the 
dichroic ARPES data and find the SAM structures. The results 
are consistent with the theoretical prediction, even in the high 
warping region. 

Bi2Te3 single crystals were grown by a vacuum Bridgman 
technique. ARPES measurement was performed at the beam 
fine 7U of UVSOR-II and HiSOR BL9A. Samples were 
cleaved in situ at 15 K and the chamber pressure was better 
than 5 x 10~^^ Torr. Photoemission spectra were acquired 
with a hemispherical photoelectron analyzer (MBS A-1 at 
UVSOR-II and Scienta R4000 at HiSOR). We used left- and 
right-circularly polarized (LCP and RCP) 8-eV photons to 
avoid signal from the bulk states. Circular polarization was 
achieved by using a MgF2 A/4 wave plate at UVSOR-II and 
by an elliptically polarized undulator at HiSOR. The energy 
resolution was about 12 meV and angular resolution was 
about 0.1°, which corresponds to a momentum resolution of 
0.0016 A-i. 

Figure 1(a) shows the experimental geometry and the 
surface Brillouin zone (BZ) of Bi2Te3. CPL comes in at 40° 
to the sample surface. Figures 1(b) and 1(c) are ARPES data 
acquired with LCP and RCP light, respectively. Dirac conelike 
dispersive bands of the TM states, especially the upper cone. 
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FIG. 1. (Color online) (a) Experimental geometry used in the study and Brillouin zone of Bi2Te3. Sample surface is parallel to the xy plane. 
Circularly polarized photons (spiral arrow) come in in the xz plane at an angle of 40° to the sample surface, (b) F-K cut taken with LCP. 
(c) Same cut taken with RCP. (d) Difference data taken with LCP and RCP. 



are observed in the data. The Dirac point is located at ~0.27 eV. 
Figure 1(d) plots the difference between LCP and RCP data. It 
shows a clear CD signal over a wide energy range. Surprisingly, 
the dichroism has very strong binding energy dependence. 
For example, in Fig. 1(d), the dichroism in the right upper 
Dirac bands switches from negative to positive as we move 
from the Dirac point toward the Fermi energy. This casts a 
contrast to the Bi2Se3 case where the dichroism has small 
energy dependence. 

One important difference between Bi2Se3 and Bi2Te3 is 
that the band warping effect is much more severe in the case 
of Bi2Te3. To see how such a warping effect plays a role in the 
dichroism, we need to inspect the CD in the momentum space 
at a constant binding energy. Figure 2(b) plots constant energy 
CD data taken with 50-meV steps between and 250 meV as 



marked by the dashed lines in Fig. 2(a). The plots show that 
the sign of the CD on a constant energy contour has a binding 
energy dependence. Figure 2(c) shows CD data at 0-, 100-, and 
250-meV binding energies. At 250 meV, which is close to the 
Dirac point, the shape of the constant energy contour is almost 
a circle. For 150-meV binding energy, the constant energy 
contour has a hexagonal shape. The Fermi contour at 0-meV 
binding energy has a strong warping effect, as manifested by 
its star shape. It is obvious that CD is also affected by the 
warping effect in Bi2Te3 and has a complex pattern. 

One can consider various possibilities as the origin of the 
dichroism. First of all, it is known that CD in ARPES can 
result from the broken chirality in the experimental geometry. 
However, in such a case, the dichroism signal is expected to 
monotonically increase as the momentum moves away from 
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FIG. 2. (Color online) (a) Illustration of Dirac cone of Bi2Te3. (b) Dichroism in the momentum space at various binding energies, (c) 
Constant dichroism map at selected binding energies of 0, 100, and 250 meV. 



the mirror plane (F-K line in this case) and cannot have the 
complex pattern shown in the leftmost panel in Fig. 2(c). 
Therefore the most probable cause for the dichroism is the 
existence of 0AM. ^"^'^^ It was previously shown that CD is 
approximately proportional to the inner product of the light 
propagation and the 0AM vectors. For TIs, Berry's phase 
refers to the relative phases of the spin-up and spin-down com- 
ponents of the wave function. In such a case, it is interesting 
to note that CD is a direct measurement of the Berry's phase 
in the k space, since the Berry's phase angle is precisely the 
angle of the spin measured from the kx axis of the Brillouin 
zone. A similar case has been reported for graphene.^^ 

CD along a constant energy contour is obtained and plotteds 
in Fig. 3. It is plotted as a function of azimuthal angle as de- 
fined in Fig. 2(c). We define the CD intensity as the normalized 
difference Icd = (LCP - RCP)/(LCP + RCP). Iqd becomes 
larger as the momentum goes away from the Dirac point and 
reaches the maximum value of about 40% at the Fermi energy. 
Such a large circular dichroism in ARPES is comparable 
to the previously reported value for Bi2Se3. Looking at the 



curves in Figs. 3(a)-3(f), it is evident that the behavior of 
Icd has very strong energy dependence. Near the Dirac point 
in Fig. 3(a), CD has an approximate form of sin^. A strong 
coupling between 0AM and SAM makes 0AM to form chiral 
states as well. In an ideal case of chiral 0AM state, the circular 
dichroism on a constant energy contour has a sin form. There- 
fore our CD data shows that 0AM has an almost ideal chiral 
structure near the Dirac point. Ideal chirality breaks down 
once the warping effect starts to come in. At 150-meV binding 
energy, which gives a hexagonal constant binding energy 
contour (weak warping effect region), Iqd is close to a sin 30 
function. By the time the energy reaches the Fermi energy, the 
warping effect is much stronger (strong warping effect region) 
and, as a result, the constant energy contour (Fermi surface in 
this case) is star-shaped. At this energy, Iqd takes a complicated 
form and requires higher order sinusoidal functions. 

As briefly mentioned above, such a complex CD pattern 
implies an equally complex 0AM structure which originates 
from the warping effect. Then what is the 0AM structure? We 
find that the CD pattern at 150 meV can be explained within 
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FIG. 3. Circular dichroism /cd as a function of at binding 
energies of (a) 250 meV, (b) 200 meV, (c) 150 meV, (d) 100 meV, (e) 
50 meV, and (f) meV. The definition for the azimuthal angle is 
shown in Fig. 2(c). 

the theoretical prediction. ^^'^^ That is, the warping effect is 
caused by higher order terms in the kp perturbation when the 
contour size becomes larger (the band has a larger k value). To 
describe the SAM structure, they discuss the in- and out-of- 
plane components separately. The in-plane component is chiral 
near the Dirac cone and turns into a wiggled chiral shape when 
the strength of the warping effect gets stronger. Meanwhile, the 
out-of-plane component is almost zero near the Dirac cone but 
becomes proportional to sin 3^ as the strength of the warping 
effect increases. 

On the other hand, circular dichroism in ARPES comes 
from the 0AM as discussed above, not from SAM. Therefore 
we need to incorporate this fact in the effective Hamiltonian 
to explain a complex CD pattern. We recently obtained 
an OAM-based effective Hamiltonian for Rashba-type band 
splitting. ^^'^^ The new Hamiltonian is given by 

HK=aK(txt)-t, (1) 

where ax is a constant that is proportional to the degree 
of asymmetric charge distribution for a given 0AM, £^ the 
surface electric field, and ~? the 0AM. The formalism for this 
new Hamiltonian simply replaces the SAM of the Rashba 
Hamiltonian by 0AM, yet it can solve the energy scale 
problem in Rashba-type split bands. Now we can write the 
kp perturbation equation for the warped Dirac cone by using 
the OAM-based terms, 

H(k) = Eo(k) + Vk(k,Ly - kyL,) + h/2(kl + kl)L, 

+ /^(^^L+-^^L_) + ... (2) 

The Hamiltonian, except that the SAM is replaced by 0AM, 
is the same as the previously reported one.^^'^^ Then we can 



expect that the complex 0AM pattern obtained from the CD 
can be explained in the same way. 

To go beyond qualitative aspect of the CD to a quantitative 
analysis, we fit the CD data in Fig. 3 and plot the result in Fig. 4. 
Data near the Dirac point plotted in Fig. 4(a) can be fitted 
with a single sin^ function, which indicates that the 0AM 
has an almost ideal helicity near the Dirac point (250 meV), 
as schematically shown in the figure. Meanwhile, fitting the 
100-meV data shown in Fig. 4(b) cannot be done with a single 
sin^ but requires an additional term of sin 3^, in the form of 
Asin ^ + Bsin 3^ . The magnitude of the sin term is somewhat 
smaller than that of the near-Dirac-point data at 250 meV. 
Here, we note that the sign of the data changes five times 
(due to the sin 3^ term), and the in-plane component alone 
cannot explain such a sign change because the chirality of the 
in-plane component cannot be reversed (that is, the in-plane 
component can change the sign change only once). Therefore, 
the sin 3^ term comes from the out-of-plane component of 
the 0AM. In fact, CPL comes in at an angle of 40° to the 
sample surface in our experimental geometry (Fig. 1) and thus 
the dichroism picks up the out-of-plane component. All things 
considered, our result is consistent with the prediction by Fu. 
In addition, it is also compatible with the recently reported 
spin-resolved ARPES data.^^'^"^ Spin-resolved ARPES data 
show that the in-plane component has a chiral structure with 
an approximately constant magnitude, while the out-of-plane 
component changes in sign and gives the sin 30 term. 

As noted above, CD at the Fermi surface has a more 
complicated shape, which implies that the 0AM texture is also 
more complex. We first note that Iqd at meV has a double- 
peak structure near = 27c/3 [arrow in Fig. 4(c)]. This feature 
has been observed in repeated measurements and thus is not 
from noise. It turns out that fitting such a structure requires 
an additional term of sin 6^. This can be explained by the 
fact that higher order terms are needed in the kp perturbation 
theory as the size of the momentum increases. In fact, the size 
of the constant energy contour of Bi2Te3 at the Fermi level is 
almost twice larger than that of Bi2Se3 . The effect from higher 
order terms is discussed in a recent theoretical work by Basak 
et al}^ on the spin texture of the strong warping case. They 
consider higher order terms to explain the nonorthogonality 
between the spin and the momentum. Nonorthogonality (5 
in the report), expressed as a function of 0, oscillates as 
sin 6^. Therefore, the sin 6^ modulation in our data may 
have come from the nonorthogonality of in-plane-direction 
0AM caused by the large Fermi surface. Besides, noting 
that the sin 60 modulation which is even function breaks the 
time-reversal symmetry, we expect that sin 60 modulation can 
be caused by extrinsic channels such as magnetic impurities or 
surface-bulk scattering. It is interesting to note that bottom 
of the bulk conduction band exists in the strong warping effect 
region. 

At this stage, it is also worth mentioning that sensitivity of 
CD to the in- and out-of-plane 0AM can be quite different. As 
discussed above, the in-plane component of 0AM at 100 meV 
has a chiral structure with an approximately constant magni- 
tude. Constant magnitude of the in-plane component implies 
that the magnitude of the out-of-plane component is small 
(otherwise, the in-plane component would be more severely 
modulated). Yet, looking at the 100-meV data in Fig. 4(b), 
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FIG. 4. (Color online) CD along the constant energy contour at (a) meV, (b) 100 meV, and (c) 250 meV. Dots and lines represent the 
data and the fit, respectively. Also drawn are schematics of the in-plane (black arrow) and out-of-plane components of the 0AM. For the 
out-of-plane components, the red segment of the contour has in-to-page 0AM while blue does out-of-page, as shown by the symbols. 



one finds that sin 3^ (which comes from the out-of-plane 
component) is as large as the sin term (in-plane component). 
The ratio between sin and sin W obtained from the fitting 
is 1:1.2. Therefore, for the photon energy we used, CD seems 
to be more sensitive to the out-of-plane component than to 
the in-plane component. We expect that such a dichotomy in 
the sensitivity has photon energy dependence because the final 
state is modified as the photon energy changes. 

In conclusion, we investigated the 0AM of Bi2Te3 by using 
ARPES with circularly polarized light. Similar to CD results 
from Bi2Se3, our data show very strong CD. We find that CD 
due to 0AM is severely affected by the warping effect in the 
band structure. The 0AM close to the Dirac point has an ideal 
chiral structure (sin^) without the out-of-plane component. 
The CD shows a sin 3^ effect in the weak warping effect 
region around 150-meV binding energy. Such 6 dependence 
is compatible with the theoretically predicted spin structure 



for the case when the out-of-plane 0AM component has a 
strong contribution to the dependence. In addition, CD gets 
an extra sin 6^ term from modulation in the in-plane 0AM 
component when the warping effect becomes very strong near 
the Fermi energy. This result is consistent with theoretical 
prediction, which considers higher order terms in our modified 
kp perturbation theory that incorporates 0AM. 
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